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ABSTRACT 

The  desorption  of  gases  from  coals  from  several 
Canadian  seams  was  studied.  Initial  rates  of  desorption 
of  carbon  dioxide,  methane  and  nitrogen  from  a  pressure 
of  50  p.s.i.g.  were  determined  for  coal  fragments  of  dif¬ 
ferent  sizes.  Total  sorption  capacity  of  the  coal  to 
carbon  dioxide,  methane  and  nitrogen  at  50  p.s.i.g.  was 
determined.  Measurements  were  made  with  the  AP  apparatus, 
and  on  the  basis  of  the  AP~  index  the  coals  were  classi 
fied  as  prone  to  outburst  or  not.  Possible  relationships 
between  total  sorption  capacity,  coal  rank  and  AP  index, 
and  also  between  the  apq-60  index,  initial  rates  of  desorp 
tion  and  coal  rank,  have  been  tested  and  analysed. 
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INTRODUCTION 


General ; -  Coal  in  an  undisturbed  seam  may  contain  a  large 
amount  of  gas  at  high  pressure.  The  emission  of  gas  from 
coal  seams  during  the  extraction  process  is  one  of  the  main 
hazards  in  coal  mining  today.  The  reduction  of  gas  pressure, 
such  as  occurs  when  the  coal  is  exposed  and  worked  at  the 
coal  face,  permits  desorption  to  take  place.  The  gas  given 
off  by  the  coal  and  associated  strata  enters  the  ventilating 
air  stream  and  pollutes  the  mine  atmosphere,  and  unless  it 
is  immediately  mixed  with  air  sufficient  to  dilute  it  to  a 
safe  level,  it  becomes  an  ignition  or  explosive  hazard. 

In  some  cases  the  desorption  of  gas  can  take  place 
very  rapidly  and  violently,  and  this  process  of  desorption 
is  described  as  a  sudden  outburst  of  coal  and  gas.  In  an 
outburst,  a  quantity  of  fine  coal  varying  from  a  few  tons  to 
several  thousand  tons  is  ejected.  Sometimes  the  fine  mater¬ 
ials  are  carried  for  several  yards,  completely  filling  the 
roadways.  In  such  a  case  ventilation  of  the  district  may  be 
cut  off  and  the  mine  atmosphere  becomes  oxygen  deficient. 
Where  such  outbursts  are  apt  to  occur  it  is  difficult  to 
continue  mining  operations,  and  production,  planning  and 
safety  of  personnel  are  greatly  affected.  The  gases  evolved 
during  outbursts  are  usually  methane,  but  in  some  cases 
carbon  dioxide  and  a  small  percentage  of  nitrogen  and  higher 
hydrocarbon  gases  are  released  (Rice,  Wilson,  Hargraves). 
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The  mechanism  of  gas  liberation  from  coal  has  attrac¬ 
ted  attention  from  the  beginning  of  this  century,  and  a  limi¬ 
ted  amount  of  research  work  has  been  done  to  establish  the 
basic  factors  which  control  the  release  of  gases  from  coal. 
Gas  pressure,  gas  concentration,  ground  pressure  and  the 
nature  of  the  coal  appear  to  be  the  main  factors  affecting 
the  liberation  of  gases  from  coal. 

The  object  of  the  work  that  is  reported  in  this  thesis 
was  to  study  certain  factors  that  are  related  to  the  release 
of  gases  from  coal  samples  from  Canadian  seams.  For  this 
purpose  studies  were  made  of  the  following: 

1.  The  initial  desorption  rate  of  gases  (CC^,  CH^,  N^) 
from  different  size  ranges  of  coal  from  various  Canadian 
sources . 

2.  The  desorption  rate  of  methane  from  coal  samples 
using  the  AP  apparatus,  which  has  been  developed  and  used  in 
Europe,  in  order  to  classify  the  different  types  of  coal 
according  to  the  AP  index. 

3.  The  total  sorptive  capacity  of  different  types  of 
coal  in  order  to  correlate  sorptive  capacity  to  the  initial 
rate  of  desorption. 

• 

Occurrence  of  Gases  in  Coal:-  It  is  considered  that  the 
gases  that  are  found  in  coal  either  have  been  formed  within 
coal  during  its  formation  or  else  have  been  introduced  from 
some  outside  source.  Their  mode  of  origin  and  their  associa¬ 
tion  with  coal  will  be  described  in  the  following  section. 
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Coal  has  been  formed  by  a  complicated  chemical  process 
from  plant  materials  and  during  this  process  some  of  the 
gases  (CH4,  CC>2 ,  higher  hydrocarbon  gases)  contained  in  the 
coal  were  formed.  The  initial  stage  in  the  formation  of  coal 
was  the  accumulation  of  vegetable  matter  which  was  then  con¬ 
verted  to  peat  under  certain  local  conditions,  then  in  se¬ 
quence  to  lignite,  and  if  maturation  continues,  to  bituminous 
coal  and  finally  to  anthracite  (Venter  and  Stassen) .  The 
basic  organic  components  of  the  vegetable  matter  are  cellu¬ 
lose,  lignin,  wax  and  resin,  in  varying  proportions  accord¬ 
ing  to  the  different  kinds  of  original  vegetable  matter. 

The  sequence  of  changes  during  coal  maturation  has 
been  described  by  several  writers  (Graham,  1949;  Venter  and 
Stassen;  Hargraves) .  Graham  (1949)  reports  that  in  the 
change  from  wood  to  lignite  with  a  carbon  content  change 
from  50  to  60%  (on  a  dry  mineral-matter  free  basis,  d.m.m.f.) 
the  main  loss  is  of  moisture,  with  release  of  a  small  amount 
of  methane  and  a  still  smaller  quantity  of  carbon  dioxide. 

In  the  conversion  of  low-rank  lignite  to  low-rank  bituminous 
coal,  i.e.,  a  change  of  carbon  from  65  to  79%,  the  chief 
loss  in  weight  is  due  to  the  loss  of  carbon  dioxide.  In  the 
next  stage  of  the  coalif ication  process,  in  which  the  carbon 
content  becomes  altered  from  79  to  86%,  he  considers  that 
there  is  a  relatively  small  production  of  methane  with  release 
of  a  larger  volume  of  carbon  dioxide.  In  the  further  change 
from  86  to  88%  carbon,  the  proportion  of  methane  formed  is 
increased,  and  during  the  change  from  88  to  92.3%  carbon  the 
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main  loss  in  weight  is  due  to  loss  of  methane.  The  final 
stage  of  change  from  semi-anthracite  to  anthracite,  92.3  to 
95%  carbon,  is  accompanied  by  a  change  in  weight  which  is 
due  almost  entirely  to  the  loss  of  methane. 

Another  source  of  carbon  dioxide  according  to  several 
writers,  including  Rice  and  Wilson,  is  from  nearby  igneous 
activity.  Carbon  dioxide  may  be  derived  from  igneous  intru¬ 
sions  and  may  have  migrated  into  the  coal  bed  through  faults 
or  joints.  Hargraves  refers  to  Russian  studies  by  Lidine, 
who  considers  the  occurrence  of  carbon  dioxide  gas  in  seams 
in  deep  levels  to  be  entirely  due  to  igneous  activity,  and 
in  the  upper  levels  to  be  due  to  permeation  of  air  down  into 
the  seam  from  the  atmosphere,  resulting  in  oxidation  of  coal 
to  produce  carbon  dioxide  and  residual  nitrogen. 

Association  of  Gas  and  Coal:-  Coal  is  a  porous  rock  and  it 
provides  a  large  internal  surface  area  for  the  adsorption  of 
gas.  According  to  several  writers,  including  Van  Krevelen 
and  Schuyer,  coal  contains  two  pore  systems,  a  macro-pore 
system  which  is  accessible  to  mercury  under  pressure,  and  a 
micro-pore  system  which,  even  at  very  high  pressures,  cannot 
be  permeated  with  mercury.  Helium  can  easily  permeate 
through  both  systems  at  room  temperature.  A  convenient  clas¬ 
sification  of  pores  according  to  their  average  width  has  been 
proposed  by  Dubinin  (1955) .  Pores  of  widths  below  20A°  are 
described  as  micropores,  those  with  width  above  200  A0  as 
macropores,  and  those  with  width  between  20  A°  and  200  A0  as 
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transitional  (or  intermediate)  pores. 

Gases  produced  during  the  coalif ication  period  may 
amount  to  several  thousands  of  cubic  feet  per  ton  (Graham, 
1949) .  The  low-rank  coals  (lignite  and  brown  coal)  which  are 
usually  found  at  shallow  depths  are  almost  always  free  from 
gases,  because  the  gases  evolved  can  usually  escape  through 
the  permeable  strata  above.  The  thicker  and  more  compact  the 
covering  deposit,  the  more  difficult  it  is  for  gases  to  es¬ 
cape.  It  is  the  normal  condition  in  deeply  bedded  coal  to 
find  hydrocarbon  gases  in  the  coal  itself  and  in  the  enclosing 
rocks.  The  strata  above  and  below  the  coal  can  also  contain 
gases  to  the  extent  of  one-twentieth  that  of  the  coal,  and 
can  act  as  a  large  reservoir  for  gas  storage  (Graham,  1922) . 

The  gases  contained  in  the  coal  can  be  held  in  the 
following  ways: 

1.  Gas  can  be  retained  in  the  micropores,  crevices, 
cracks  and  bedding  planes  of  the  coal  in  more  or  less  a  free 
gaseous  state,  so  that  it  may  be  given  off  rapidly  when  these 
openings  are  intersected  by  the  mine  workings. 

2 .  Gas  can  be  adsorbed  in  large  quantities  on  the  inter¬ 
nal  surface  of  coal  which  is  provided  by  the  microporous 
nature  of  coal,.  The  process  of  adsorption  is  different  than 
that  of  absorption  of  gases.  When  a  gas  is  allowed  to  come 
to  equilibrium  with  a  solid  or  liquid  surface,  the  concentra¬ 
tion  of  gas  molecules  is  always  greater  in  the  immediate 
vicinity  of  the  surface  than  in  the  free  gas  phase,  regard¬ 
less  of  the  nature  of  the  gas  or  surface.  The  process  by 
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which  this  surface  excess  is  formed  is  called  adsorption. 
According  to  Young  and  Crowell,  the  phenomenon  of  adsorption 
is  like  a  condensation  of  vapour  to  form  a  film  of  liquid. 

It  is  to  be  distinguished  from  absorption,  which  involves 
bulk  penetration  of  the  gas  into  the  structure  of  the  solid 
by  some  process  of  diffusion.  In  the  case  of  coal  and  gas 
it  is  difficult  to  distinguish  between  those  processes,  so 
the  term  'sorption'  is  generally  used,  which  embraces  both 
types  of  phenomenon. 

The  adsorbed  gas  forms  a  very  thin  layer  on  the  exter¬ 
nal  and  internal  surfaces  of  coal.  In  addition,  some  gases 
will  be  held  by  mechanical  compression  in  the  micropore 
spaces.  Both  forms  of  gas  are  liberated  when  the  ambient 
pressure  is  reduced,  especially  when  the  coal  is  broken  up 
by  crushing  and  grinding. 
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HISTORICAL  SURVEY  OF  PREVIOUS  INVESTIGATIONS 


The  problem  of  gas  release  in  the  working  of  coal 
faces  had  been  realised  before  the  beginning  of  this  century. 
Much  of  the  work  done  on  this  problem  has  been  in  an  effort 
to  obtain  further  knowledge  of  outbursts  and  their  incidence. 

A  limited  amount  of  work  has  been  carried  out  in  the  labora¬ 
tory  to  investigate  the  coal-gas  relationships  and  to  study 
the  various  factors  which  are  responsible  for  the  desorption 
of  gases  from  coal.  A  short  review  of  certain  selected 
aspects  of  previous  work  is  given  in  the  following  section, 
including  sorption  studies  and  the  effect  of  pressure  and 
temperature,  rate  of  desdrption  studies,  measurement  of  the 
quantity  and  pressure  of  gas  in  situ,  and  measurement  of 
permeability. 

Sorption  Studies;-  The  phenomenon  of  sorption  under  differ¬ 
ent  conditions  has  been  studied  by  several  workers.  Graham 
(1917)  studied  the  solubility  of  different  gases  in  coal  and 
he  showed  that  although  coal  is  a  solid  substance,  it  can 
dissolve  considerable  quantities  of  gases  and  in  general, 
the  more  easily  liquifiable  the  gas  the  more  readily  will  it 
be  dissolved  by  coal.  He  found  that  the  sorptive  capacity 
of  coal  for  carbon  dioxide  was  greater  than  for  methane . 

Later,  in  1922,  he  reported  studies  on  the  sorption  of  various 
gases  in  coal  at  pressures  up  to  30  atmospheres. 
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Venter  and  Stassen  refer  to  studies  made  by  Coppen, 
in  which  Coppen  showed  that  sorption  of  gas  in  coal  at  first 
increases  in  proportion  to  pressure,  but  at  high  pressure 
the  gas  is  sorbed  to  a  much  more  restricted  degree. 

Moffat  and  Weale  studied  the  sorptive  capacity  of 
several  types  of  British  coal  to  methane  at  pressures  up  to 
1,000  atmospheres.  Moffat  and  Weale  measured  total  sorption 
(methane  sorbed  by  the  lump  volume  of  coal)  and  apparent  sorp¬ 
tion  (methane  sorbed  only  on  the  internal  surface) .  In 
their  studies  they  found  that  total  sorption  usually  in¬ 
creases  directly  with  pressure  at  first  and  then  more  slowly 
above  100  -  150  atmospheres.  The  highest  value  of  total 
sorption  varied  according  to  the  type  of  coal  from  28  to  70 
ml./gm.  (1,000  to  2,500  cu.  ft. /ton).  Similarly,  apparent 
sorption  also  increases  directly  with  pressure  to  100  -  150 
atmospheres,  where  total  sorption  in  this  case  reaches  a 
maximum  of  20  to  61  ml./gm.  and  of  varying  sharpness,  accord¬ 
ing  to  the  coal.  Moffat  and  Weale  also  observed  the  relation¬ 
ship  of  sorptive  capacity  to  coal  rank  at  different  pressures 
and  noted  that,  at  any  particular  high  pressure,  the  sorptive 
values  for  coals  of  carbon  content  between  85  and  92%  were 
lower  than  those  for  coals  having  higher  or  lower  carbon  con¬ 
tent. 

Gunther  in  reporting  some  of  the  studies  made  at 
CERCHAR  (Centre  d' Etudes  et  Recherches  des  Charbonnages  de 
France)  noted  that: 

1.  The  sorption  of  CH4,  CC>2  and  N2  by  coal  is  reversible 
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and  reproducible; 

2.  The  sorption  isotherm  has  the  approximate  shape  of  a 
regular  hyperbola,  and  that; 

3.  CC>2  is  sorbed  more  than  CH^  which  in  turn  is  sorbed 
more  than  N2 ,  but  this  difference  is  less  noticeable  at  high 
pressures . 

In  the  previous  paragraphs  it  has  been  noted  that  the 
quantity  sorbed  in  general  increases  in  proportion  to  in¬ 
creasing  pressure  up  to  a  certain  pressure,  and  that  the  in¬ 
crease  in  sorption  is  not  so  great  at  high  pressures. 

The  effect  of  temperature  and  pressure  on  sorption  of 
gases  has  been  studied  by  several  workers.  Graham  (1917) 
studied  the  effect  of  temperature  on  the  sorption  of  various 
gases  and  found  that  sorption  decreases  with  rise  of  tempera¬ 
ture.  Coppen  in  his  studies  also  found  a  decrease  in  sorption 
with  a  rise  of  temperature.  Gunther  reported  that  an  increase 
of  1°  C.  in  temperature  reduces  the  sorption  capacity  of  coals 
by  about  0.8%. 

Rate  of' Desorption  Studies:-  Several  workers  in  the  past 
have  studied  the  rate  of  desorption  by  different  means.  Some 
have  studied  the  rate  of  desorption  from  the  initial  period 
of  desorption  until  there  was  no  further  yield  of  gas,  while 
others  have  studied  the  rate  of  desorption  only  in  the  first 
seventy  second  period.  Botham  and  Patching  (1959)  studied 
the  rate  of  emission  of  gas  from  coal  quickly  taken  from  bore¬ 
holes,  using  a  portable  apparatus  and  procedure  devised  by 
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Hargraves  of  the  University  of  Sydney,  Australia.  Fragments 
of  fresh  coal  from  a  borehole  were  quickly  placed  in  tubes 
that  were  then  connected  to  manometers ,  and  the  rates  at 
which  gas  was  emitted  were  found  by  observing  the  rates  of 
fluid  displacement o  Initial  rates  of  gas  desorption  were 
observed  to  be  much  greater  from  samples  from  outbursting 
seams.  Later  Patching  (1966)  made  preliminary  studies  of 
the  rates  of  desorption  under  controlled  conditions,  especi¬ 
ally  during  the  initial  few  seconds  after  release  of  pressure. 

Belin  refers  to  the  index  method  which  was  perfec¬ 
ted  by  Somnier  from  the  above  method  introduced  by  Hargraves. 

A  sample  of  10  gm.  of  coal  (quickly  taken  from  a  borehole  at 
the  coal  face  at  a  depth  of  2  -  3  meters)  is  screened  to  the 
range  of  0.5  -  0.8  mm.  and  is  quickly  poured  into  a  sample 
holder  and  connected  to  a  manometer  where  the  gas  emitted 
displaces  a  fluid  at  atmospheric  pressure.  The  index 
represents  the  volume  of  gas  released  by  the  sample  from  the 
35th  to  the  70th  second  after  sampling. 

In  Australia,  Bolt  and  Innes  studied  the  rate  of  emis¬ 
sion  of  gases  from  types  of  coal  in  the  laboratory  and  sug¬ 
gested  that  diffusivity  was  not  a  dominant  factor  in  the 
mechanism  of  outbursts.  Hargraves  reported  studies  on  the 
rate  of  emission  of  gases  from  different  types  of  coal  under 
different  conditions,  and  he  observed  that  for  coals  of  equal 
total  sorption  capacity,  but  unequal  permeability,  the  rate- 
of-emission  curve  for  the  more  permeable  coal  has  a  steeper 
initial  gradient  and  a  more  rapidly  flattening  gradient  than 
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non-permeable  coal. 

In  France,  Belgium,  Poland  and  the  U.S.S.R.,  rates  of 
desorption  studies  have  been  done  in  detail  on  the  basis  of 
the  index  and  on  what  is  called  the  aP  index.  This  AP 
index  is  obtained  using  equipment  and  procedures  which  have 
been  slightly  modified  from  the  original  method  of  Ettinger 
of  the  U.S.S.R.  A  quantity  of  coal  is  degassed  for  1^  hours 
and  then  saturated  for  1%  hours  with  methane  at  a  pressure 
of  one  atmosphere.  It  is  then  connected  to  a  receptacle  — 
initially  evacuated  —  of  known  cubic  content.  As  the  gas 
is  released,  the  pressure  in  the  receptacle  gradually  rises, 
and  the  increase  in  pressure  (AP)  after  a  certain  time  is 
recorded.  The  index  AP  defines  the  initial  speed  of  desorp¬ 
tion.  A  more  detailed  account  of  this  method  will  be  given 
in  the  experimental  section  and  in  the  appendix. 

Gunther  (1965)  mentions  that  a  relationship  between 
the  AP  index  and  the  index  in  the  French  mines  has  been 
established  as 


C 


AP 

870 


where  C  is  gas  concentration  in  the  coal  in  equilibrium  with 
pressure  P  (gas  pressure  in  the  fissures) . 


Measurement  of  Quantity  and  Pressure  of  Gases  In  Situ:- 

Coal  seams  which  were  previously  worked  by  hand  are  now 
usually  machine  cut.  With  high  rates  of  advance  of  coal 
faces  the  pressure  gradient  of  the  gas  in  the  coal  and  sur¬ 
rounding  strata  tends  to  be  greater,  giving,  consequently, 
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an  increased  rate  of  evolution  of  gas.  Knowledge  of  the 
content  and  pressure  of  gas  in  the  coal  would  be  helpful  in 
assessing  the  possible  dangers  from  emission  of  gases.  A 
number  of  attempts  have  been  made  in  the  past  to  measure  the 
content  and  pressure  of  gas  in  situ,  Graham  (1938)  measured 
the  residual  gas  content  from  a  fresh  lump  of  coal  (1^  x  1 
x  1  in.  taken  quickly  from  the  face  and  put  into  a  container) 
and  the  total  evolution  of  firedamp  was  found  to  be  500  cc . 
for  22.7  gm.  of  anthracite.  He  also  measured  the  gas  pres¬ 
sure  in  situ  and  recorded  seam  pressures  up  to  40  atmos¬ 
pheres.  Botham  and  Patching  (1959)  measured  the  gas  content 
ahead  of  the  coal  face  in  the  Upper  Marsh  seam  (Canmore 
Mines  Ltd.)  and  they  found  that  there  was  an  increase  in  gas 
content  of  the  coal  with  distance  into  the  solid.  Some  gas 
(50  to  115  cu.  ft./ton  at  N.T.P.)  was  found  in  the  coal  on 
the  exposed  faces,  and  the  amount  generally  increased  to 
about  400  -  450  cu.  ft. /ton  at  depths  of  12  to  20  feet. 

Cervik  has  reported  gas  pressures  in  coal  seams  of  up  to  35 
atmospheres  in  American  mines. 

Measurement  of  Permeability:-  Although  there  are  undoubt¬ 
edly  many  factors  that  affect  the  emission  of  gas  from  coal, 
permeability  is  one  of  the  most  important.  The  permeability 
of  coal  has  been  measured  by  several  workers.  Graham  (1920) 
measured  the  permeability  of  coal  to  various  gases  and  he 
reported  that  it  was  impossible  for  gas  to  pass  through  coal 
if  the  coal  were  free  from  cracks  and  breaks.  Recent  studies 


13 


of  permeability  measurements  include  those  of  Botham  (1958) , 
Skochinsky  and  Khodot,  Haung  and  Shelton,  Patching  (1965) 
and  Gunther.  Botham,  like  Graham,  found  that  there  was  no 
apparent  flow  of  gas  through  solid  unfissured  coal  specimens. 
Patching  (1965)  summarised  some  of  his  permeability  measure¬ 
ments  on  coals  that  were  deformed  or  disturbed  to  various 
degrees : 

1.  Permeability  depends  essentially  on  flow  through 
fractures  or  microcracks. 

2.  Appreciable  variations  in  permeability  can  be  expec¬ 
ted  between  samples  from  the  same  seam. 

3.  Permeability  is  drastically  reduced  by  the  applica¬ 
tion  of  confining  pressure. 

4.  Gas  permeability  is  reduced  as  mean  gas  pressure 
increases,  all  other  factors  being  constant. 

5.  Permeability  is  not  noticeably  affected  by  gas  sorp¬ 
tion,  but  does  depend  on  the  size  of  the  gas  molecules. 


THEORY 


Discussion  of  the  Various  Proposed  Theoretical  Equations  of 

Diffusion; -  The  emission  of  gases  from  coal  depends  on 
several  factors  including  pressure  and  concentration  of  gas, 
permeability  of  coal,  fracturing,  coal  particle  size,  and 
the  diffusion  coefficient.  On  the  basis  of  Darcy's  and 
Fick's  laws,  several  different  theories  have  been  developed 
to  explain  the  mechanism  of  gas  emission  from  coal.  A  short 
review  of  the  proposed  theories  and  their  applications  will 


be  discussed  in  the  following  section. 

In  1856,  Henry  Darcy  formulated  a  law  for  fluid  flow 
in  porous  media  which  may  be  expressed  as  follows; 


q' 
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where  q'  =  volume  flow  rate  (in  cm. /sec.) 
k  =  permeability  (in  darcys) 
y  =  gas  viscosity  (in  centipoises) 

A  =  cross  sectional  area  (in  cm?) 

P  =  pressure  (in  atmospheres) 

1  =  length  (in  cm.) 


(1) 


The  driving  force  for  this  mode  of  transport  is  the  pressure 


gradient.  The  movement  of  gas  through  the  fracture  or  pore 
system  can  be  considered  to  take  place  according  to  the  above 
law  in  those  conditions  where  gas  can  easily  flow  through 


the  fractured  spacing  of  the  coal  bed. 

Fick,  using  Fourier's  heat  conduction  equation  as  an 
analogy,  postulated  a  mathematical  theory  of  diffusion  in  an 
isotropic  substance  on  the  hypothesis  that  the  rate  of 
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transfer  of  diffusing  substance  through  unit  area  of  a 
section  is  proportional  to  the  concentration  gradient  meas¬ 
ured  normal  to  the  section,  ice., 


'  p,  i  _  b  dC 

q  "  dl 

where  q*  =  rate  of  transfer  per  unit  area 
D  =  diffusion  coefficient 
C  =  concentration  of  gas  in  solid  coal 
1  =  length 


(2) 


The  driving  force  for  this  type  of  transport  is  a  concentra¬ 
tion  gradient .  From  equation  (2)  we  can  derive  the  differ¬ 
ential  equation 

2 

dc  _  D  d  C  , 

dt  “  775 

dl 


Expressions  (2)  and  (3)  are  usually  referred  to  as  Pick's 
first  and  second  laws  of  diffusion. 

Lipayev  (1956,  equation  33)  proposed  an  expression 
for  the  release  of  gas  from  a  sphere  of  coal,  on  the  basis 
of  conservation  of  mass  of  gas,  and  he  assumed  that  flow 
followed  the  linear  law  of  Darcy „  Conservation  of  mass  re¬ 
quires  that  the  mass  flowing  into  the  volume  element  minus 
the  mass  flowing  out  of  the  element  in  a  given  time  interval 
is  equal  to  accumulation .  He  took  an  ideal  case,  and  assumed 
that  the  initial  piece  of  coal  and  all  its  parts  resulting 
from  crushing  had  a  spherical  form.  His  analysis  applied  to 
a  case  in  which  a  piece  of  coal  was  suddenly  ejected  from  a 
seam,  in  which  the  initial  concentration  of  the  gas  was 
assumed  to  be  Po  atmospheres . 

The  rate  of  gas  liberation  in  this  case  was  given, 
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according  to  Lipayev,  by  the  equation 

2 

dQ  4 ^  _  k  ,  —  7t  at,  ... 

=  RT  Po  Ri  “  9  exp  (— — 2—)  (4) 

Ri 

where  ^  -  volume  flow  rate 

Po  =  initial  gas  pressure 

R]_  =  radius  of  the  spherical  particle 

k  =  coefficient  of  gas  permeability  of  coal 

y  =  viscosity  of  methane 

g  =  gas  density 

a  =  constant 

t  =  time 

R  =  gas  constant 
T  =  temperature  °k 

For  a  particular  sample  of  coal  with  a  given  initial  gas  con¬ 
tent,  and  with  a  given  size  of  particle,  this  equation  may 
be  written  simply  as 


H  =  ke'Xt  (4a) 

and  integrating  this  equation  for  the  total  yield  of  gas 
at  time  t 

Qt  =  Q„  (1  -  e_At)  (5) 

Equation  (5)  satisfies  the  boundary  conditions  needed  for 
diffusion  of  gas  from  coal,  i.e., 

when  t  =  0  Q  =  0 

and  when  t  =  00  Q  =  finite  =  Q 

oo 

Lipayev  suggested  that  the  equation  might  be  checked 
experimentally  by  measuring  the  expenditure  of  methane  from 
a  conic  sample  when  a  large  gauge  tap  is  suddenly  opened. 
Considering  the  cone  as  an  element  of  a  sphere,  it  should  be 
possible  to  obtain  by  calculation  the  rate  of  gas  escape 
from  a  sphere  of  coal.  Tests  could  be  made  on  spherical 
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fragments  of  coal.  However,  experimental  verification  of 
this  equation  is  not  available  in  the  literature. 

Patching  (1959)  tried  to  fit  this  equation  with  his 
experimental  results  (which  were  done  with  fragments  of  coal 
and  did  not  include  data  for  the  first  minute)  but  the  above 
equation  did  not  fit  the  observed  data. 

Bolt  and  Innes  used  the  following  equation  to  fit 
their  experimental  studies  of  gas  emission  from  coal  parti¬ 
cles  . 

Qt  =  (1  -  Ae“At)  (6) 


where  Qt  =  volume  of  gas  desorbed  at  any  time  t 

=  ultimate  gas  yield 

t  =  time 

A,X  =  constants 

This  equation  is  an  abbreviated  form  of  the  solution  provi¬ 
ded  by  Crank  of  Fick's  second  law.  The  solution  of 
equation  (3)  is  given  by  Crank  on  the  assumption  that  the 
diffusion  coefficient  (D)  is  always  constant  and  that  the 
particles  are  of  spherical  size  of  radius  R-^ .  According  to 
Crank  (equation  6.20  p.  86)  the  total  amount  of  diffusing 
substances  that  has  entered  or  left  a  sphere  at  any  time  t 
is  given  by 
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(7) 


In  order  to  determine  the  value  of  D  the  following 

approximation  may  be  used.  From  the  first  three  terms  of 

t  Qt 

equation  (7) ,  the  value  of  —  0  for  =  0,5  is  given  by 
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and  with  approximation  of  this  equation,  it  is  possible  to 
calculate  D  from  the  relation 

72 

D  =  0.035  (8a) 

r0.5 

where  tA  c  is  the  time  at  which  =  0.5 

0°5  Q~ 

Bolt  and  Innes  tested  equation  (6)  against  their 
experimental  data  and  they  concluded  that  this  equation  is 
valid  after  the  first  minute  of  desorption,  but  noted  that 
in  the  first  minute  the  actual  rate  of  desorption  was  higher 
than  given  by  the  equation.  It  is  apparent  that  this  equa¬ 
tion  does  not  satisfy  the  initial  boundary  condition  unless 
A  =  1,  whereas  they  found  values  of  A  to  be  about  0.6.  With 
the  help  of  equation  (8a)  they  calculated  diffusion  coeffi¬ 
cients  for  several  types  of  coal. 

Other  workers,  including  Nelson  and  Walker,  and  Nandi 
and  Walker,  have  made  use  of  equation  (7)  for  their  diffu¬ 
sion  studies  of  gas  from  coal,  and  have  calculated  various 
parameters  including  the  diffusion  coefficient  and  the  acti¬ 
vation  energy.  Many  investigators  in  the  past  have  apparent¬ 
ly  assumed  the  diffusion  coefficient  to  be  constant  but 
recently  some  modifications  have  been  suggested  by  Russian 
workers.  Khodot,  Yanovskaya  and  Premyslev  refer  to 
a  modification  of  equation  (7) ,  as  proposed  by  Yanovskaya 
and  Timofeev,  namely 
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Qt  =  Q< 
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where  r  = 


IT 


Def  ^0.5 


Rl2  A 


Def  =  effective  diffusion  coefficient 
A  =  coefficient  of  adsorption  of  methane  in  coal 

=  (Vn  is  the  content  of  free  gas  in  the  pores) 
v0  u 

Yanovskaya  and  Timofeev  found  experimentally  that  the 
diffusion  coefficient  is  not  constant  irrespective  of  gas 
pressure  and  coal's  fissuration,  but  it  is  dependent  on 
these  two  parameters.  Their  equation,  like  that  of  Bolt  and 
Innes,  does  not  fit  the  initial  boundary  conditions. 

Because  equation  (7)  does  not  converge  rapidly  for 
small  values  of  time,  an  alternate  solution  of  Pick's  second 
law  equation  for  small  values  of  time  has  been  derived 
(Crank,  equation  6.22,  p.  87)  in  the  different  form  of 


Qt  ,  ,Dt  -  h  ,  -H  .  ,  ”  •  -  nRl,  ,Dt 
TT-  =  6  ( — T)  [tt  +  2  E,  lerfc  ]  -  3 — . 

Q~  ^  "-1  /Dt  Rx' 


(10) 


By  an  approximation  of  equation  (10)  Barrer  and  Brook 
provided  a  different  solution  for  the  second  law  of  Pick  as 
follows : 

Q, 


t  _  2s  /Dt 

Q  v  t/  tt 
*^00 

where  s  =  total  external  surface  of  the  grains 
v  =  total  volume  of  the  grains 


(11) 


Equation  (11)  has  been  used  at  CERCHAR  for  diffusion 
studies.  It  is  reported  by  Gunther  that  when  the  dimensions 
of  a  piece  of  coal  exceed  the  distance  between  fissures,  the 
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relationship  above  is  no  longer  valid.  With  the  use  of 

equation  (11)  and  measurements  made  on  fragments  of  5  -  10  y 

size  for  which  s  and  v  were  determined  by  direct  examination 

with  a  microscope,  the  diffusion  coefficient  was  found  to  be 

-10  2 

of  the  order  of  10  cm./sec.  for  coal  with  CH^  at  room 
temperature.  If  the  grain  size  is  so  large  that  the  influ¬ 
ence  of  fissures  is  noticeable,  then  he  noted  that  use  of 
equation  (11)  yields  a  higher  value  of  diffusion  coefficient 
than  if  tests  are  made  on  small  granules  of  the  same  coal. 
Gunther  also  reports  that  if  measurements  of  s  and  v  are 
done  by  the  BET  method,  supposing  that  nitrogen  is  not 
sorbed  at  all  at  77  °K  (which  is  also  debatable)  in  this 
case  results  lead  to  undervaluation  of  D. 

As  a  second  type  of  approach,  an  equation  has  been 
derived  at  CERCHAR  on  the  basis  of  Darcy's  law,  assuming 
that  coal  is  a  fissured  medium  in  which  gas  moves  only  in 
the  fissures  according  to  Darcy's  law.  The  final  equation 
is  reported  by  Gunther  as  follows: 


kP  dP  , 
div  —  j—  grade  = 
y  dc 


°°  2 

/  (t-DR(t)dr 

0  3t2 


(12) 


where  k  =  permeability  of  fissured  material  in  darcies 
P  =  gas  pressure  in  the  fissure 

c  =  gas  concentration  in  the  coal  in  equilibrium  with 
pressure  P,  expressed  in  cm^  (adjusted  at  N.T.P.) 
per  cm3  of  material 

r(t)  =  the  function  of  retardment.  This  function 
increases  with  time,  being  0  for  t  -  0  and 
1  for  t  “ 


According  to  Gunther,  at  the  present  time  there  is  no 
particular  solution  for  the  above  equation,  but  it  may  be 
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solved  someday  by  numerical  methods.  Experimental  verifi¬ 
cation  of  equation  (12)  is  not  available  in  the  literature 
but  Gunther  reports  that  some  practical  solutions  have  been 
provided  at  CERCHAR  in  several  particular  cases. 

To  the  present,  there  appears  to  be  no  single  equation 
which  can  adequately  describe  the  mechanism  of  gas  emission 
from  coal.  It  is  proposed  to  use  several  of  the  above  equa¬ 
tions  to  interpret  experimental  data  of  the  release  of  gas 
from  coal,  and  the  application  of  these  equations  will  be 
described  in  the  discussion  section  of  this  thesis. 
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EXPERIMENTAL  STUDIES 


General : -  Three  different  types  of  tests  were  carried  out 
on  coal  from  several  different  sources.  Studies  were  made 
of  the  following: 

1.  The  initial  desorption  rate  of  gases  (C02,  CH N2) 
from  different  size  ranges  of  coal  from  various  Canadian 
sources . 

2.  The  desorption  rate  of  methane  from  coal  samples 
using  the  AP  apparatus  in  order  to  classify  the  different 
types  of  coal  according  to  this  index. 

3.  The  total  sorptive  capacity  of  the  same  samples  of 
coal  in  order  to  correlate  sorptive  capacity  to  the  initial 
rate  of  desorption. 

The  source  of  coal,  the  description  of  the  apparatus, 
and  the  procedures  employed  are  given  in  the  following. 

The  Source  of  Coal  Samples:-  For  the  studies  described  in 
this  thesis  samples  from  three  different  mining  companies 
were  used.  A  visit  was  made  in  May  1967  to  the  Canmore  Mine 
Limited,  and  lump  coal  samples  were  obtained  from  17  Slope 
(at  the  side  of  an  outburst  zone)  and  from  the  face  of  41 
Slope  of  the  Upper  Marsh  seam  in  No.  3  Mine,  and  from  the 
Wilson  seam.  Later,  in  September  1967,  another  visit  was 
made  to  the  Canmore  Mines  Limited  and  coal  samples  were 

l 

taken  from  channels  (from  roof  to  floor)  from  3  Crosscut, 
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37  Slope,  and  from  3  Crosscut,  41  slope,  of  the  Upper  Marsh 
seam.  Figures  1,  la  and  lb  show  the  locations  where  samples 
were  taken  from  the  Upper  Marsh  seam. 

At  Canmore,  there  have  been  outbursts  in  the  Upper 
Marsh  seam,  but  not  in  the  Wilson  seam.  "The  coal  at  Can- 
more  is  of  Lower  Cretaceous  age,  and  the  coal  bearing  form¬ 
ation  occurs  as  an  assymetric  syncline  with  the  southwestern 
limb  partly  overturned  by  the  thrust  of  the  older  Paleozoic 
formations  from  the  southwest.  Tectonic  movements  during  the 
folding  of  the  coal  basin  caused  many  small  folds  and  faults 
in  the  coal  seams  and  enclosing  strata.  The  thickness  of 
the  Upper  Marsh  seam  varies  from  about  2  to  3*5  meters,  be¬ 
cause  of  these  folds  and  faults  and  also  because  of  the 
original  irregular  deposition  of  the  coal  deposits." 
(Patching,  1966)  . 

Coal  (samples  1  and  lc)  from  the  Upper  Marsh  seam 
(Canmore  Mines)  is  very  soft  and  friable.  It  is  hard  to 
distinguish  between  bright  (vitrain,  clairain)  and  dull 
(durain,  fusain)  layers.  The  coal  of  sample  2  (Wilson  seam) 
is  tighter  than  that  from  the  Upper  Marsh  seam  and  alternate 
bright  and  dull  bands  are  visible. 

Coal  was  also  provided  from  different  seams  by  Crows- 
nest  Industries  Limited  and  McIntyre  Mines  Limited. 

"The  coal  seams  in  the  Crowsnest  area  were  laid  down 
in  the  Lower  Cretaceous  period,  and  the  coal  bearing  forma¬ 
tion  in  the  Western  district  occurs  in  a  basin  with  the 
center  covered  by  a  great  thickness  of  younger  rocks  except 
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FIGURE  1A  LOCATION  OF  SAMPLES 

UPPER  MARSH  SEAM,  CANMORE 
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FIGURE  IB  LOCATIONS  OF  SAMPLES 

UPPER  MARSH  SEAM,  CANMORE 
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where  cut  by  sharp  valleys.  There  has  been  a  variable 
amount  of  faulting  and  folding  which  has  caused  the  coal  to 
be  often  sheared  and  very  friable."  (Patching).  In  the 
past,  outbursts  have  occurred  in  the  Crowsnest  area,  although 
not  in  the  seams  which  are  being  mined. 

Coals  from  the  Crowsnest  Industries  look  alike, 
except  samples  3d  and  3e.  In  most  of  the  samples  alternate 
bands  of  bright  and  dull  coals  are  visible.  Not  much  mineral 
matter  (fusain)  is  present.  In  general,  the  coals  from  this 
district  have  a  tighter  structure  than  coal  from  the  Upper 
Marsh  seam  or  the  Wilson  seam.  However,  sample  3e  (C  seam) 
is  soft  and  friable,  and  has  a  less  tight  structure  than  the 
rest  of  the  samples  from  Crowsnest  Industries.  Bands  are 
visible  in  the  coal  from  this  seam,  and  there  are  more 
durain  and  fusain  parts  of  coal  than  vi train  and  clairain. 
Coal  samples  from  C  seam  and  those  from  the  Upper  Marsh 
seam  look  much  alike  except  that  the  coal  from  C  seam  has  a 
somewhat  tighter  structure  than  the  coal  from  the  Upper 
Marsh  seam. 

Little  information  is  available  to  the  author  about 
the  coal  deposits  of  McIntyre  Mines  Limited.  All  of  the 
samples  (4a,  4b  and  4c)  from  this  mine  have  similar  appear¬ 
ance.  Bands  are  visible  and  some  patches  of  fusain  are 
present.  The  coal  samples  of  McIntyre  Mines  in  general  are 
very  hard  and  have  tighter  structure  than  the  samples  from 
Canmore  or  the  Crowsnest  district. 


Table  1  shows  the  sources  of  the  coal  used  for  the 
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TABLE  1:  PROXIMATE  ANALYSIS  OF  THE  COAL  SAMPLES 
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Analysis  of  the  coal  samples  was  done  by  the  Research  Council  of  Alberta 
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experimental  studies  and  the  coal  analysis  of  several  samples 
is  also  given. 

Sample  Preparation:-  Lumps  of  coal  as  received  from  the  mines 
were  crushed  to  -V  and  screened  to  the  following  sizes 
(Tyler  mesh):  -3+4,  -4+8,  -8+14,  -14+28,  -28+48,  -48+65, 
-65+100  and  -100+200. 

For  measurements  with  the  AP  apparatus  only  samples 
of  0.25  -  0.50  mm.  (-32+60  mesh)  were  used. 
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TEST  PROCEDURE  NO.  1 

This  series  of  experiments  was  carried  out  to  deter¬ 
mine  the  initial  rates  of  desorption  of  gas  from  samples  of 
coal  which  had  been  charged  with  gas  at  50  p.s.i.g.  The 
apparatus  included  a  vacuum  pump,  a  source  of  gas  (C02,  CH^, 
N2  and  He)  ,  a  manifold  to  which  six  metal  sample  holders 
could  be  connected,  an  expansion  chamber  (29  cm.  long  and 
8.97  cm.  internal  diameter),  a  differential  pressure  trans¬ 
ducer  (0  -  1  p.s.i.  range  with  10  volt  D.C.  excitation)  and 
a  chart  recorder.  The  sample  holders  could  contain  6.5  gm. 
of  coal  fragments  of  any  size.  They  consisted  of  cups 
threaded  at  the  top  and  provided  with  detachable  caps  which 
could  be  screwed  tightly  against  an  O-ring  seal.  Each  cap 
was  provided  with  two  on/off  valves. 

The  equipment  is  shown  schematically  in  Figure  2 . 

Operational  Procedure;-  Portions  of  coal  of  a  given  size 
(e.g.,  -3+4  mesh)  and  source,  weighing  6.5  gm. ,  were  placed 
in  the  cups.  They  were  evacuated  for  45  minutes  and  then 
charged  with  gas  (CC>2,  CH^ ,  N2  or  He)  at  a  pressure  of  50 
p.s.i.g.  for  twenty  minutes  to  one  hour.  This  period  was 
found  to  be  enough  for  the  sorption  of  gas  by  the  coal 
samples.  The  supply  of  gas  was  then  closed  off  and  the 
exhaust  valve  of  the  sample  holder  suddenly  opened,  allowing 

i 

the  rapid  release  of  gas  into  the  expansion  chamber.  The 
increase  of  pressure  in  the  expansion  chamber  due  to  the 
addition  of  the  desorbed  gas  was  detected  by  the  pressure 
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FIGURE  2  DIAGRAM  OF  DESORPTION  APPARATUS 


■ 


33 


transducer  and  recorded  by  the  chart  recorder.  This  pres¬ 
sure  build-up  was  recorded  only  for  the  initial  two  minutes. 
Two  or  three  such  traces  of  the  pressure  build-up  were  ob¬ 
tained  for  each  size  and  sample  of  coal  with  the  same  gas. 
Similar  determinations  were  made  for  each  size  of  each  sample 
of  coal  using  methane  and  nitrogen. 

Corrections  were  applied  for  the  free  volume  in  each 
cup,  including  the  intergranular  and  free  pore  space  that 
was  occupied  by  gas.  For  this,  samples  were  evacuated,  char¬ 
ged  with  helium  at  a  pressure  of  50  p.s.i.g.  for  ten  minutes 
and  then  discharged,  and  the  traces  were  recorded  in  the 
same  way  as  before.  Since  the  helium  was  not  sorbed  by  the 
coal,  the  trace  obtained  was  due  to  the  helium  in  the  free 
spaces,  and  this  provided  a  quantity  that  could  be  subtrac¬ 
ted  from  the  traces  provided  by  the  other  gases. 

In  this  experiment  the  curve  produced  by  the  chart 
recorder  indicated  the  quantity  of  gas  desorbed  at  any 
given  time,  and  by  suitable  calibration  the  actual  quantity 
of  gas  could  be  calculated.  Comparison  of  emission  rates 
could  be  directly  made  by  comparison  of  plots. 
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TEST  PROCEDURE  NO,  2  (AP  APPARATUS) 

This  apparatus  (Figure  3)  was  made  to  the  same  dimen¬ 
sions  as  the  Ap  apparatus  designed  by  CERCHAR.  There  are 
eight  metal  cups  of  duraluminum  which  can  each  hold  three 
grams  of  coal.  Each  cup  is  provided  with  a  groove  for  an 
O-ring  on  the  top  so  the  flange  of  a  glass  valve  can  make  a 
good  seal  with  the  cup.  A  ball  and  socket  joint  clip  is 
used  to  keep  each  cup  and  valve  together.  By  turning  the 
knob  of  the  stopcock,  the  interior  volume  of  any  cup  can  be 
connected  either  to  a  front  or  rear  manfold.  Samples  are 
charged  by  helium  or  methane  through  the  front  manifold. 

The  exhaust  (rear)  manifold  can  be  connected  to  a  vacuum 
pump,  and  a  manometer  is  also  connected  to  the  rear  manifold. 
The  manometer  is  set  in  front  of  a  millimetric  graduation  of 
25  cm.  in  length.  The  total  volume  of  one  cup  plus  rear 

manifold  (with  the  mercury  at  the  same  level  in  both  limbs 

3  3 

of  the  U-tube)  is  26  cm. ,  of  which  7  cm.  is  for  each  cup  and 

3 

19  cm.  for  the  exhaust  manifold. 

Operational  Procedure:-  Three  grams  of  coal  of  -32+60  mesh 
were  placed  in  each  cup.  Samples  were  degassed  for  1%  hours 
and  then  charged  with  helium  at  a  pressure  of  one  atmosphere 
(760  mm.)  Each  cup  in  turn  was  then  connected  with  the  rear 
manifold  in  which  a  vacuum  was  re-established  between  each 
measurement.  The  difference  in  the  manometer  levels  were 
noted  for  each  sample  in  mm.  of  Hg  (called  P^) .  The  value 
of  P^  is  a  measure  of  the  contribution  of  free  gas  which  is 
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FIGURE  3:  Photograph  of  the  AP  apparatus. 
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not  sorbed  by  the  sample  but  simply  held  under  atmospheric 
pressure  in  the  dead  volume  formed  by  the  volume  of  the  coal 
pores  and  the  intergranular  spaces. 

The  helium  was  now  evacuated  from  all  samples  and 
lines,  and  then  methane  was  introduced  by  the  front  manifold 
at  atmospheric  pressure  for  a  period  of  1^  hours.  This  period 
is  considered  to  be  enough  to  saturate  the  samples  with 
methane . 

After  saturation,  the  vacuum  in  the  exhaust  manifold 
was  re-established  and  for  each  sample  the  indices  V  ^  and 
were  noted  as  follows:  a  cup  to  be  tested  was  connected  to 
the  exhaust  manifold  for  10  seconds  and  then  the  tap  was 
closed,  and  the  difference  in  the  levels  of  manometer  was 
noted  as  mm.  of  Hg.  Again  at  the  end  of  35  seconds,  i.e. 
45  seconds  after  the  beginning  of  the  test,  the  tap  was 
again  connected  to  the  rear  manifold  for  15  seconds  and  then 
the  tap  was  closed,  and  again  the  difference  in  manometer 
levels  was  noted  as  P^  mm.  of  Hg,  thus  corresponding  to  a 
desorption  period  of  60  seconds.  Two  or  three  portions  of 
each  type  of  coal  were  tested  and  then  characterised  by  the 
AP  indices. 


APq_£ o  values  were  calculated  as  follows: 
AP0-10  =  P2  “  P1 


AP10-60  P3  P2 

AP0-60  =  AP0-10  +  AP10-60  P3  "  P1 


Measurements  with  the  AP  apparatus  can  be  done  to  an 


accuracy  of  ±4  mm.,  i.e.,  the  AP  indices  determined  for  a 
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given  sample  of  coal  usually  agree  within  four  units. 

For  further  details  of  the  equipment  and  procedure 
involved,  refer  to  Appendix  II. 
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TEST  PROCEDURE  NO.  3 

With  the  apparatus  described  in  the  first  series  of 
tests,  it  was  only  practical  to  measure  the  volume  of  gas 
desorbed  up  to  2  minutes  but  not  to  measure  the  total  amount 
of  gas  desorbed  from  a  coal  specimen.  To  obtain  more  com¬ 
plete  information  about  the  desorption  mechanism,  it  was 
necessary  to  know  the  ultimate  quantity  of  gas  that  would 
be  desorbed.  Therefore  a  third  series  of  experiments  was 
carried  out  to  determine  the  total  sorption  capacity  of  the 
different  coals. 

Two  similar  sets  of  equipment  were  used.  Each  set 
included  a  source  of  gas  (C02,  CH^ ,  N2  and  He),  a  metal  cell 
(used  as  in  Experiment  1),  a  100  cc .  burette  with  three-way 
stopcock,  and  a  levelling  bulb. 

The  equipment  is  shown  schematically  in  Figure  4 . 

Operational  Procedure:-  Portions  of  -100+200  mesh  and  of 
-14+28  mesh  fragments  from  the  same  sample  of  coal,  each 
weighing  6.5  gm. ,  were  placed  in  different  cells.  The  sam¬ 
ples  were  evacuated  for  45  minutes  and  then  charged  with 
carbon  dioxide  at  a  pressure  of  50  p.s.i.g.  at  22°  C.  for 
2  to  18  hours.  The  supply  valve  was  then  closed  and  the 
exhaust  valve  opened  so  that  the  desorbing  gas  could  go 
into  the  burette  and  displace  the  mercury.  The  mercury 
levels  in  the  burette  and  the  levelling  bulb  were  adjusted 
to  the  zero  mark  of  the  burette  at  the  beginning,  and  as  gas 
flowed  into  the  burette  the  levelling  bulb  was  always  kept 
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FIGURE  4  DIAGRAM  OF  DESORPTION  APPARATUS 
TO  MEASURE  SORPTION  CAPACITY 
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in  pace  with  the  mercury  level  of  the  burette,  and  the  cor¬ 
responding  volume  and  time  readings  were  noted  until  there 
was  no  further  desorption  of  gas.  The  burette  readings  could 
be  estimated  to  the  nearest  0.2  ml.  for  the  0  -  100  ml. 
burette.  Similar  determinations  were  made  for  the  same  sizes 
and  samples  of  coal  using  methane  and  nitrogen. 

For  each  different  type  of  coal,  the  free  volume 
occupied  by  the  gas  in  the  cell  and  in  the  pore  spaces  was 
determined  by  charging  the  samples  with  helium  gas  at  50 
p.s.i.g.  for  10  minutes  and  by  measuring  the  volume  of 
helium  released  on  expansion „  This  quantity  was  subtracted 
from  the  corresponding  quantity  of  other  gases  used  in  the 
test  in  order  to  determine  the  volume  of  gas  that  was  desor¬ 


bed. 
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RESULTS  AND  DISCUSSION 


Total  Sorption  Capacity:-  The  total  amount  of  gas  desorbed 
from  a  sample  of  coal  in  the  present  studies  is  not  the 
amount  of  gas  that  was  sorbed  by  that  coal,  because  the 
samples  at  the  end  of  the  test  yre re  still  connected  to  a 
supply  of  the  gas  at  atmospheric  pressure  and  therefore  some 
residual  gas  undoubtedly  remained  in  the  samples.  Hence, 
results  are  not  directly  comparable  with  those  of  other  in¬ 
vestigators  who  measured  the  amounts  sorbed.  However,  it 
is  considered  that  the  test  conditions  of  desorption  repre¬ 
sent  in  a  reasonable  way  the  conditions  that  do  exist  at  the 
actual  face  of  an  advancing  mine  heading. 

Table  2  shows  the  sorption  capacities  of  various  sam¬ 
ples.  The  values  shown  are  the  average  of  two  to  three  tests 
and  the  average  deviation  from  the  mean  value  was  found  to 
be  +5  cu.  ft. /ton.  It  can  be  seen  from  Table  2  that  the 
total  desorbed  volume  appears  to  differ  for  different  part¬ 
icle  sizes  of  the  same  coal  sample.  In  order  to  determine 
if  the  coal  in  the  -14+28  fragments  had  the  same  sorption 
capacity  as  the  coal  that  broke  to  -100+200  fragments,  sev¬ 
eral  samples  of  -14+28  mesh  size  fragments  were  crushed  to 
-100+200  mesh  and  the  gas  sorption  re-determined.  The 
results  obtained  from  these  determinations  are  shown  with 
asterisks  in  Table  2  and  they  are  approximately  equal  to 
previous  values  of  sorption  capacities  obtained  from  -100+200 
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mesh.  It  appears  that  the  capacity  of  the  coal  was  the  same 
for  both  sizes,  but  that  the  large  sizes  obviously  were  not 
fully  saturated. 

It  was  found  that  there  was  virtually  no  significant 
difference  in  the  amount  of  gas  desorbed  from  the  -100+200 
mesh  fragments,  whether  one  hour,  two  hours  or  eighteen 
hours  was  allowed  for  the  saturation  period.  Therefore  it 
can  be  concluded  that  the  -100+200  mesh  sizes  were  saturated 
at  test  conditions;  however,  the  time  allowed  was  not  enough 
for  saturation  of  the  -14+28  mesh  sizes,  although  in  many 
cases  these  samples  were  charged  with  gas  for  a  period  of 
ten  to  eighteen  hours. 

Table  2  shows  that  the  total  sorption  capacity  differs 
to  some  extent  between  the  different  samples  of  coal.  Fig¬ 
ure  5  was  plotted  to  determine  if  there  exists  any  relation¬ 
ship  between  the  total  sorption  capacity  and  the  coal  rank. 

It  is  not  possible  to  draw  any  conclusion  from  this  plot 
regarding  the  variation  of  sorption  capcity  with  coal  rank. 

Tests  were  not  made  to  determine  the  difference  be¬ 
tween  total  sorption  and  apparent  sorption  at  the  test  pres¬ 
sure  of  50  p.s.i.g.  Recently  it  was  determined  by  Jolly, 
Morris  and  Hinsley  that  there  was  no  appreciable  difference 
between  the  total  sorption  and  the  apparent  sorption  at  pres¬ 
sures  of  less  than  10  atmospheres.  Therefore  in  the  test 
studies,  no  distinction  is  made  between  the  two  types  of 
sorptions . 

In  the  present  study  the  least  and  the  highest  values 
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of  total  desorption  quantity  were  observed  for  different 
samples  having  87%  carbon  content  (d.nwm.f.)  Moffat  and 
Weaie  found  a  minimum  sorption  by  coal  having  carbon  con¬ 
tents  between  85  and  92% 0  Jolly,  Morris  and  Hinsley  obser¬ 
ved  a  minimum  of  sorption  between  85  and  90%  carbon  content. 
They  also  suggested  that  it  is  unlikely  that  a  general  cor¬ 
relation  of  methane  sorption  with  other  parameters,  such  as 
rank,  can  be  obtained  until  the  sorption  capacities  of  a  very 
large  number  of  coals  are  obtained. 

The  total  desorption  of  carbon  dioxide  was  found  to 
be  greater  than  that  of  methane  and  the  amount  of  methane 
in  turn  was  greater  than  that  of  nitrogen.  The  ratio  of  the 
total  amount  of  desorbed  carbon  dioxide  to  methane  was  found 
to  average  1.7:1  and  the  range  was  from  lt4  to  2.4  and  the 
ratio  of  carbon  dioxide  to  nitrogen  averaged  3=8:1.  This 
carbon  dioxide  to  methane  ratio  is  in  approximate  agreement 
with  the  results  of  Ettinger,  Chaplinsky,  Lamba  and  Adamov, 
who  found  in  their  studies  that  the  ratio  of  carbon  dioxide 
to  methane  was  from  1.5  -  1.7  to  1. 

They  noted  that  the  difference  in  the  sorption  capa¬ 
city  of  coals  corresponds  to  the  higher  temperature  of  liqui- 
f action  of  carbon  dioxide  (-78.5°  C.)  as  compared  with  me¬ 
thane  (-161.5°  C.)  and  nitrogen  (-195°  C.)  Generally  it 
has  been  observed  that  at  any  given  temperature,  a  certain 
sorbent  can  sorb  more  readily  those  gaseous  substances  that 
condense  at  high  temperature  (Maron  and  Prutton) . 

In  this  series  of  experiments  the  effects  of  different 
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pressures,  temperatures  and  moisture  contents  on  the  total 
sorption  capcity  were  not  studied.  All  samples  where  char¬ 
ged  with  gas  at  50  p.Soi.g.,  and  the  ambient  temperature  was 
22°  C.  ±  0.5°  Co  during  the  tests.  It  has  been  shown  by 
several  investigators,  including  Moffat  and  Weale,  and  Jolly, 
Morris  and  Hinsley,  that  the  sorption  capacity  of  coal  to 
gas  increases  nearly  linearly  with  pressures  up  to  100  atmos¬ 
pheres.  Therefore,  knowing  the  sorption  capacity  at  a  given 
pressure,  a  rough  estimate  of  gas  sorption  can  be  made  at 
other  pressures. 

The  effect  of  temperature  has  been  studied  by  several 
workers.  In  general  they  have  found  that  the  sorption  capa¬ 
city  decreases  with  rise  of  temperature.  Recently,  Jolly, 
Morris  and  Hinsley  have  found  that  a  rise  of  10°  C.  in  temp¬ 
erature  decreases  the  sorption  capacity  approximately  by  20%. 
These  other  studies  indicate  the  temperature  variations  that 
occurred  during  the  tests  that  are  herein  reported  were  not 
large  enough  to  cause  any  significant  variations  in  the  test 
data. 

Rates  of  Desorption  Studies:-  The  rates  of  desorption  of 
gases  (C02,  CH^  and  N2)  from  different  sizes  and  samples  of 
coal  are  shown  in  Figures  6  to  16.  These  curves  are  obtained 
by  averaging  two  or  three  duplicate  runs  and  tracing  the 
charts  at  actual  scale.  In  each  case  there  is  a  definite 
pattern  of  desorption.  Table  3(a)  shows  the  quantities  of 
gases  (C02 ,  CH^  and  N2)  desorbed  at  27  seconds.  Here  again 
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TABLE  3(a):  QUANTITY  OF  GAS  DESORBED  AT  27  SECONDS  FROM 
INITIAL  PRESSURE  50  p.s.i„g.  (-100+200  mesh  size) 
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it  is  seen  from  this  table  that  the  amount  of  carbon  dioxide 
desorbed  at  27  seconds  is  greater  than  that  of  methane  and 
the  amount  of  methane  in  turn  is  greater  than  that  of  nitro¬ 
gen  at  the  same  time.  The  same  theoretical  explanation 
applies  as  discussed  for  the  total  sorption  of  the  various 
gases.  It  is  worth  noting  that  the  ratio  of  C02  to  CH^ 
emitted  at  27  seconds  is  on  the  average  3.50:1  in  comparison 
to  1.7:1  at  infinite  time  (i.e.,  when  desorption  is  complete). 
This  higher  ratio  of  carbon  dioxide  to  methane  in  the  begin¬ 
ning  of  the  desorption  is  likely  due  to  the  higher  sorption 
capacity  of  the  coal  for  carbon  dioxide  than  for  methane  at 
any  given  pressure.  Therefore  a  higher  concentration  of 
carbon  dioxide  will  be  present  and  released  in  the  early 
stages  of  desorption  as  compared  to  the  quantity  of  methane 
that  would  be  released. 

Values  of  the  diffusion  coefficient  were  calculated 
by  using  equation  (8a) 


D  «  0.035  Tr± — 

0.5 

for  eleven  samples  of  -100+200  mesh  size,  and  are  shown  in 

Table  3(b).  The  three  samples  from  McIntyre  Mines  Limited, 

having  carbon  content  of  77%  (d.m.m.f.  basis),  had  the  low- 

-7  2 

est  values  of  D  (0.20  x  10  cm  )  while  samples  from  Canmore 

sec 

Mines  Limited,  having  carbon  content  of  88%  had  high  values 
-7  2 

(0.92  x  10  cm  )  of  D.  Sample  3e  (C  seam,  Crowsnest)  with 
sec 

a  carbon  content  of  67%  has  also  the  same  high  value  (0.92  x 
-7  2 

10  cm  )  of  D.  It  can  be  seen  that  there  is  a  wide  range  of 
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variation  of  D  with  coal  rank,  and  the  variation  of  D  does 
not  follow  a  definite  trend  with  the  increasing  or  decreasing 
carbon  content.  Therefore  it  is  not  possible  to  draw  any 
relationship  between  the  diffusion  coefficient  and  the  coal 
rank,  but  in  general  the  samples  (e.g.  from  McIntyre)  which 
are  tightest  in  appearance  have  shown  the  lowest  values  of 
the  diffusion  coefficient.  The  values  obtained  for  D  are  in 
general  agreement  with  those  of  Bolt  and  Innes. 

From  Figures  6  to  15  it  can  be  seen  that  the  initial 
rate  of  desorption  is  much  higher  from  samples  composed  of 
fine  granules.  This  follows  from  the  fact  that  for  samples 
of  a  given  weight  of  fine  size  particles,  there  is  a  larger 
total  surface  area  available  than  for  samples  of  coarse 
grain  particles  of  the  same  total  weight.  Another  explana¬ 
tion  is  arrived  at  from  equation  (6) : 

Qt  =  a.  d  -  Ae"xt> 

2 

where  X  =  tt  D 


It  follows  from  the  above  equation  that  for  a  constant  value 
of  Qot  and  A  and  t,  Qt  will  decrease  with  increasing  particle 
size  (R^) .  It  is  possible,  using  the  above  relationship,  to 
calculate  the  quantity  of  gas  (CC^)  desorbed  at  any  parti¬ 
cular  time  (e.g.  27  seconds)  from  -14+28  and  -100+200  mesh 
sizes  and  then  to  compare  this  quantity  with  experimental 
results.  This  has  been  done  for  two  samples,  as  shown  in 
the  following  table. 
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Comparison  of  Computed  and  Experimental  Desorbed 


Volume 

!  of  Carbon 

Dioxide  at 

27  seconds 

• 

Sample 

Number 

-14+28 

-100+200 

Computed 
values 
cu. ft/ ton 

Expt' 1 . 
values 
cu. ft/ ton 

Computed 
values 
cu. ft/ton 

Expt'  1 . 
values 
cu. ft/ton 

1  c 

112 

117 

227 

237 

2 

116 

122 

211 

220 

However,  from  Figure  16  it  is  obvious  that  there  is 
not  much  difference  in  the  rates  of  desorption  from  different 
size  particles  of  this  particular  sample.  The  coal  (sample 
1,  Upper  Marsh  seam)  is  from  an  outbursting  seam  and  it  is 
so  fractured  and  sheared  that  there  is  not  much  difference 
in  the  effective  fragmentation  despite  differences  in 
particle  size,  and  therefore  rates  of  desorption  are  approx¬ 
imately  the  same  from  different  size  range  particles. 

Figures  17  and  18  are  comparative  graphs  of  the  rates 
of  emission  of  carbon  dioxide  and  methane  for  one  size  range 
(-100+200)  from  different  samples  of  coal.  It  is  evident 
that  the  initial  rate  of  desorption  is  in  many  instances 
appreciably  different  for  different  samples  of  coal  of  the 
same  size.  Samples  1  and  lc  (Upper  Marsh  seam) ,  2  (Wilson 
seam)  and  3d  and  3e  (Crowsnest)  show  high  rates  of  desorp¬ 
tion,  whereas  samples  4a,  4b  and  4c  (McIntyre)  show  the 
slowest  rates  of  desorption.  From  these  observations  it  is 
evident  that  certain  samples  are  more  permeable  than  others, 
where  there  is  easier  escape  for  gas  molecules  from  the  in¬ 
terior  of  the  coal. 
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TABLE  3  (b) :  CALCULATED  VALUES  OF  D  AND  THE  COMPUTED  VALUES 
OF  A  AND  A  (Eqn.  6)  AND  OF  X  (Eqn.  5)  FOR  ELEVEN  SAMPLES 
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TABLE  4:  COMPUTED  VALUES  OF  Qt  (for  C02)  FROM  EQUATIONS  (5)  AND  (6) 
_ SAMPLE  #1,  41  SLOPE ,  LfPPER  MARSH  SEAM,  CANMORE _ 

Particle  Sizes 

-100+200  mesh  -14+28  mesh 
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Attempts  were  made  to  fit  several  of  the  theoretical 
equations  to  the  data  from  test  runs  on  the  different  sam¬ 
ples  of  coals.  In  equation  (5)  there  is  one  unknown.  A,  and 
there  are  two  unknowns  in  equation  (6),  A  and  A.  Qot  (total 
sorption  capacity)  was  determined  experimentally  for  all 
samples  and  these  values  of  have  been  shown  in  Table  2. 

A  separate  computer  program  was  run  for  each  equation  to  get 
the  best  possible  values  of  A  (equation  5)  and  of  A  and  A 
(equation  6)  on  the  basis  of  the  Least  Squares  method.  Val¬ 
ues  of  Qt  were  computed  from  the  resulting  equations  in 
order  to  be  compared  with  the  experimentally  determined  val¬ 
ues  of  Qt«  The  computer  programs  are  shown  in  Appendix  3a 
and  3b.  Computed  value  of  A  and  of  A  and  A  for  eleven  sam¬ 
ples  are  shown  in  Table  3(b).  The  experimental  data  and  the 
corresponding  computed  values  of  Qt  from  equation  (5)  and 
(6)  for  sample  1  (41  slope,  Upper  Marsh  seam)  for  the  -14+28 
and  -100+200  mesh  sizes  are  shown  for  different  values  of  time 
in  Table  4.  A  graph  of  Qt  (experimental)  versus  log  t  for 
-100+200  mesh  was  plotted  for  each  of  the  other  samples 
(Figures  19  to  29) .  The  corresponding  computed  values  for 
equations  (5)  and  (6)  were  also  plotted  on  the  first  three 
graphs.  It  is  obvious  from  Figures  19,  20  and  21  that  equa¬ 
tion  (5)  does  not  fit  the  experimental  data,  but  equation 
(6)  does  fit  the  experimental  data  within  reasonable  limits. 
Therefore,  the  computed  values  obtained  from  equation  (6) 
only  were  plotted  for  the  rest  of  the  samples  (Figures  22 
to  29).  It  can  be  seen  that  for  small  values  of  time,  i.e.. 
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for  the  initial  period  of  desorption,  there  are  considerable 
deviations  between  the  experimental  and  the  computed  values 
even  with  equation  (6)  .  This  demonstrates  that  equation  (6) 
does  not  satisfy  the  entire  range  of  the  desorption  period. 

It  is  worth  noting  from  Figures  19  to  29  that  for  the 
initial  period  of  desorption  there  is  apparently  a  straight 
line  relationship  between  Q  (experimental)  and  log  t.  How¬ 
ever,  the  total  amount  of  gas  that  is  emitted  from  a  given 
portion  of  coal  must  be  a  finite  quantity,  while  a  straight 
line  plot  on  a  semi-log  scale  has  no  finite  limit  and  more¬ 
over,  it  has  no  beginning.  Thus  there  appears  to  be  no 
theoretical  justification  for  such  a  relationship  and  the 
apparently  straight  line  relationship  must  be  only  approxi¬ 
mately  true  over  a  portion  of  the  desorption  period. 

Equations  (5)  and  (6)  were  also  tested  by  comparing 
computed  values  of  Qt  against  the  volumes  of  carbon  dioxide 
desorbed  from  samples  of  different  size  fragments  during  the 
same  time  interval  (27  seconds) .  A  tabulation  of  the  experi¬ 
mental  and  computed  values  of  Qt  from  equations  (5)  and  (6) 
for  sample  1  (Upper  Marsh  seam)  is  shown  in  Table  5.  (Simi¬ 
lar  tabulations  were  prepared  for  the  other  samples,  but  are 
not  included  in  this  thesis.)  Graphs  of  Qt  versus  log 
(size)  for  several  samples  of  coal  were  plotted  (Figures  30 
to  32) .  Computed  values  of  obtained  from  equations  (5) 
and  (6)  were  also  plotted  on  Figure  30  (a)  and  (b) .  Here 
again  computed  values  from  equation  (5)  showed  great  devia¬ 
tions  from  the  experimental  values.  From  Figures  30  and  31 
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it  can  be  seen  that  the  computed  values  obtained  from  equa¬ 
tion  (6)  are  very  close  to  the  experimental  data  for  several 
samples.  However,  in  some  cases  this  equation  does  not  fit 
the  whole  range  of  experimental  data,  as  for  sample  3c  (C 
seam,  Figure  31b) ,  and  sample  2  (Wilson  seam,  Figure  32) . 

It  is  also  obvious  from  Figures  30,  31  and  32  that  the  com¬ 
puted  values  obtained  from  equation  (6)  show  some  deviation 
from  the  experimental  data  for  particle  sizes  greater  than 
892  microns.  It  can  be  seen  from  Figure  32  that  the  yield 
of  gas  for  sample  1  (41  slope,  Upper  Marsh  seam)  is  almost 
independent  of  size  greater  than  200  microns.  This  is  out¬ 
burst  coal  and  highly  fractured.  In  the  light  of  the  above 
tests,  equation  (6)  appears  to  be  valid  for  a  certain  range 
of  particles  and  after  a  certain  period  of  diffusion. 
Equation  (7)  as  it  stands, 

2  2 

_  _  rn  6  00  1  ,-n  it  Dtv  i 

Q^.  -  [1  -  — 2  e  exP  ( - 2 — )  1  (7) 

tt  n=l  n 

was  not  tested  against  the  experimental  data.  It  was  stated 
by  Bolt  and  Innes  that  to  a  first  approximation  equation  (7) 
may  be  written  as 

2 

TT  Dt 

Qt  =  Q„  U  -  0.6  e  ~n 

and  this  approximation  is  equation  (6)  as  they  used  it. 

-7  2 

They  calculated  that  for  t  =  1  minute  and  D  =  10  cm. /sec. 
the  second  term  in  the  summation  in  equation  (7)  is  only 
one-tenth  of  the  first  term.  It  has  already  been  indicated 
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in  the  previous  discussions  that  equation  (6)  does  not  des¬ 
cribe  the  initial  period  of  desorption.  When  t  is  small  the 
second  and  third  term  etc.  of  equation  (7)  are  not  negligible 
and  the  approximation  for  equation  (6)  is  therefore  not 
valid.  However,  no  direct  computer  programming  was  done  to 
attempt  to  fit  the  unabbreviated  equation  (7)  to  the  experi¬ 
mental  values. 

Equation  (11) ,  as  developed  by  Barrer  and  Brook, 

^t  _  2s_  /Dt 

Q  V  /  7T 

00 

was  also  tested  against  the  experimental  data.  A  plot  of  Q 

w 

versus  /t  was  drawn  for  several  samples  (Figures  33  to  37) . 

It  can  be  seen  that  there  is  an  approximate  straight  line 
relationship  between  Qt  and  /t  for  the  initial  period  of 
time,  and  so  the  experimental  data  agrees  reasonably  with 
the  results  to  be  expected  from  this  equation.  However, 
for  large  values  of  t,  this  relationship  does  not  hold  true 
and  the  equation  can  only  be  considered  valid  for  the  initial 
period  of  time.  This  is  evident  from  the  nature  of  the  equa¬ 
tion  because  the  total  sorption  capacity  of  any  coal  is  a 
finite  quantity  at  infinite  time  whereas  from  equation  (11) 
a  finite  value  of  Q  is  not  reached.  It  is  also  evident 

^oo 

from  Figures  33  to  37  that  for  t  =  0,  Qt  on  extrapolation 
on  a  straight  line  will  be  greater  than  zero,  and  this  pro¬ 
jection  is  not  in  agreement  with  equation  (11)  nor  with  the 
initial  conditions  of  the  tests.  Therefore  it  can  be  con¬ 
cluded  that  equation  (11)  satisfied  the  experimental  data 
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in  the  initial  period  of  desorption  except  for  the  inter¬ 
cept  on  the  ordinate  when  t  =  0  but  is  not  valid  for  large 
values  of  time. 

Recently  it  was  proposed  by  Mullin  and  Berkowitz  that 
isothermal  rates  of  disengagement  of  CO  and  CH^  from  anthra¬ 
cite  bituminous  and  sub-bituminous  coal  chars  at  temperatures 
between  570  and  670°  C.  follow  a  pseudo-first  order  kinetics 
law  but  the  rates  of  disengagement  indicate  that  the  initial 
gas  discharge  is  governed  by  two  concurrent  control  mechan¬ 
isms.  Upon  plotting  Q  against  /t  as  expected  from  equation 
(11) ,  they  found  that  the  data  appeared  to  plot  along  two 
successive  straight  lines  intersecting  at  some  point,  thereby 
suggesting  that  two  stages  of  diffusion  were  involved.  The 
tests  of  gas  desorption  carried  out  for  the  present  thesis 
are  not  directly  comparable  to  those  of  Mullin  and  Berkowitz 
because  their  studies  were  based  at  high  temperatures  where 
formation  of  CO  and  CH^  is  possible  within  the  decomposing 
mass.  However,  from  Figures  33  to  37,  there  is  a  similar 
pattern  of  two  straight  lines  indicating  two  stages  of  dif¬ 
fusion,  but  it  should  be  noted  here  that  there  is  no  possi¬ 
bility  of  the  formation  of  CO  and  CH^  at  the  ambient  tempera¬ 
ture  of  22°  C. 

It  is  therefore,  concluded  that  none  of  the  equations 
fit  the  whole  range  of  the  experimental  data  of  desorption 
of  gas  from  coal  fragments,  and  it  is  possible  that  the  actual 
mechanism  of  desorption  is  more  involved  than  that  on  which 
these  equations  are  based. 
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AP  Measurement:-  Table  6  shows  the  APrt  values  determined 
-  0-6  0 

for  different  samples  of  coal.  Replicate  portions  of  twelve 
samples  were  tested  by  CERCHAR  as  a  check,  and  their  results 
are  also  shown  in  Table  6.  Several  graphs  were  plotted  of 
the  APq_60  index  against  coal  rank,  total  desorption  capa¬ 
city  and  initial  rate  of  desorption  (Figures  38,  39  and  40) . 
From  Figure  38  it  is  not  possible  to  draw  any  relationship 
between  the  APQ_g0  index  and  coal  rank.  This  agrees  with 
other  workers  including  Belin  and  Vandeloise,  who  have  de¬ 
termined  that  the  APq_^q  index  is  highly  dependent  on  the 
fissuration  of  coal  but  not  on  other  factors  as  total  sorp¬ 
tion  or  coal  rank.  From  Figure  39  there  is  a  slight  indi¬ 
cation  of  increase  in  APq_^q  values  with  increasing  total 
sorption  capacity.  It  is  also  evident  from  Figure  40  that 
the  APq_60  values  increase  with  high  rates  of  initial  gas 
desorption.  In  general,  samples  1,  lc  and  2  from  Canmore 
showed  high  rates  of  desorption  (100  -  108  cu.  ft. /ton  of 
methane  in  a  minute)  and  the  corresponding  APo-60  were 
between  70  and  93.  Samples  3d  and  3e  from  Crowsnest  also 
showed  high  rates  of  desorption  and  the  corresponding  apq_6q 
values  were  higher  than  the  indices  of  the  rest  of  the 
samples  from  the  Crowsnest.  Samples  4a,  4b  and  4c  from 
McIntyre  showed  slow  rates  of  desorption  and  the  correspon¬ 
ding  APq-60  values  were  obtained  in  the  range  of  11  -  15. 


Interpretation  of  the  AP^_^  Values:-  Vandeloise 

reports  the  classification  of  AP  index  values,  applicable  to 
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TABLE  6:  OBSERVED  APQ_60VALUES  FOR  SEVERAL  COAL  SAMPLES 


Coal 

Sample 

Type  of 

Location 

AP0-60  (U- 

of  A. ) 

AP0 

-60 

Number 

LOal 

Reading 

Average 

(CERCHAR) 

1 

Upper  Marsh 
41  slope 
face 

41  slope 
face 

90,92,95,90 

93 

87, 

87 

1A  (i) 

37  slope 

3  xcut 
(Upper  Marsh) 

0"  -  15" 
from  roof 

75,77,77 

76 

72, 

72 

lA(ii) 

II 

15"  -  36" 
from  roof 

43,41,45 

43 

1 A  ( i  i  i ) 

II 

36"  -  56" 
from  roof 

35,34,39 

37 

33, 

32 

lA(iv) 

II 

56"  -  72" 
from  roof 

55,55,55,50 

53 

lA(v) 

II 

72"  -  88" 
from  roof 

66,64,64,68 

66 

0"  -  13" 

lA(vii) 

II 

from  roof 
east  end 

96,94,92,94 

94 

96, 

97 

of  face 

0"  -  15" 

lA ( viii) 

II 

from  roof 
east  end 

74,74,76 

75 

of  face 

IB  ( i ) 

41  slope 

3  xcut 

0"  -  15" 
from  roof 

76,78,74 

76 

67, 

68 

lB(ii) 

II 

15"  -  36" 
from  roof 

56,56 

56 

IB (iii) 

II 

36"  -  52" 
from  roof 

41,40,39,39 

40 

36, 

30 

IB  (iv) 

II 

52"  -  72" 
from  roof 

41,42,44 

43 

1C 

17  slope 
(Upper  Marsh) 

17  slope 
side  of 
outburst 

71,71,69 

70 

2 

Wilson  seam 

— 

70,69,71,69 

70 

62, 

60 

3  (a) 

# 1  seam 
south  mine 

- 

21,21,24,21 

22 

13, 

18 

3  (b) 

Balmer 
#1  mine 

- 

17,22,21 

20 

xz 
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TABLE  6  (continued) 


Coal 

Sample 

Type  of 

Location 

AP0-60  <U- 

of  A. ) 

AP0-60 

Number 

LOal 

Reading 

Average 

(CERCHAR) 

3  (c) 

#7  seam 
strip  mine 

- 

17,23 

20 

3  (d) 

Balmer 

north 

- 

48,50,46,46 

48 

50,  51 

3  (e) 

C  seam 

- 

55,50,53,55 

53 

49,  54 

4  (a) 

#1  mine 
(top) 

- 

14,14,16,14 

15 

13,  14 

4  (b) 

#1  mine 
(bottom) 

- 

14,14,12 

13 

4  (c) 

#2  mine 

— 

10,12,12,12 

11 

10,  14 
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Belgian  coals  with  sudden  outbursts  of  methane,  as  follows: 


1 

° 

<3 

60 

Degree  of  Danger 

0  - 

15 

Coal  not  prone  to  outburst 

15  - 

30 

Slightly  suspected  coal 

30  - 

45 

Suspected  coal 

45  - 

60 

Dangerous  coal 

> 

60 

Highly  dangerous  coal 

The  above  classification  is  based  on  the  measurements 
of  a  sufficient  number  of  samples  (at  least  100) ,  well  dis¬ 
tributed  over  the  working  place,  and  according  to  such  meas¬ 
urements  the  susceptibility  of  the  panel  to  outburst  may  be 
fairly  assessed. 

It  can  be  seen  from  Table  6  that  the  APq_^q  values 
vary  in  a  range  from  10  to  95 0  Only  the  last  three  samples 
fall  in  the  safe  category,  while  the  remaining  samples  are 
in  the  ranges  of  slightly  suspected  coal  and  highly  danger¬ 
ous  coal.  Sample  2  (Wilson  seam) ,  according  to  the  apq_60 
classification,  lies  in  the  category  of  highly  dangerous 
coal,  but  according  to  our  present  knowledge  there  have  been 
no  sudden  outbursts  from  this  particular  seam.  There  are 
probably  several  reasons  for  this.  Firstly,  there  is  not 
much  pressure  due  to  overburden  on  this  seam,  as  mining 
there  is  done  under  less  than  700  feet  of  cover  at  present. 
Secondly,  this  coal  seam  is  not  as  faulted  and  fissured  as 
the  Upper  Marsh  seam.  In  confirmation  of  this,  it  is  sug¬ 
gested  by  several  writers  including  Vandeloise  that  the 

APrt  alone  cannot  be  used  to  characterise  the  risk  com- 
0-6  0 

pletely.  It  is  also  necessary  to  take  into  account  the 
amount  of  gas  present  in  the  seam  and  the  state  of  stress 
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in  the  coal.  Another  point  should  be  taken  into  considera¬ 
tion,  namely  that  the  experimental  Apq_60  values  as  deter¬ 
mined  are  not  the  average  values  of  general  areas  but  are 
only  the  values  of  samples  from  a  particular  place  in  a  seam. 
Samples  3d  and  3e  (Crowsnest)  also  come  in  the  category  of 
dangerous  coal,  but  again  there  is  no  report  of  outbursts 
from  these  particular  seams. 

Samples  1,  1A  (i) ,  1A  (viii)  and  1C,  all  from  Upper 
Marsh  seam,  show  very  high  values  of  the  APo_go  ^ndex'  an& 
indicate  highly  dangerous  coal.  Outbursts  have  been  repor¬ 
ted  from  time  to  time  from  the  Upper  Marsh  seam  where  these 
samples  were  taken.  It  is  worth  noting  that  samples  of  the 
1A  and  IB  series  have  widely  differing  values  of  the  APo_go 
index,  being  greatest  at  the  upper  section  (0"  -  15")  of  the 
seam  and  smaller  for  the  rest  of  the  section,  although  all 
of  these  samples  have  nearly  the  same  sorption  capacity  at 
50  p.s.i.g.  and  also  the  same  coal  rank.  This  agrees  with 
the  visual  inspection  of  the  seam  in  which  it  was  observed 
that  the  upper  section  of  the  seam  at  the  sampling  locations 
was  more  sheared  than  the  rest  of  the  sections.  This  agrees 
also  with  the  observation  of  Belin  and  Loison  (1964)  and 
others  that  the  AP  index  often  varies  appreciably  between 
different  parts  of  the  same  seam. 

It  can  be  surmised  on  the  basis  of  the  above  obser¬ 
vations  that  there  was  a  greater  possibility  of  outbursts 
occurring  in  the  upper  sections  of  the  seam  where  the  highest 
values  of  the  APn  index  were  observed  than  in  the  lower 
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sections  of  the  seam  at  these  sampling  locations.  This  cor¬ 
responds  with  general  experience  in  the  mine  that  the  great¬ 
est  displacement  of  coal  during  outbursts  has  occurred  near 
the  roof  of  the  seam. 


■ 

\ 

' 


CONCLUSIONS 


From  the  experimental  work  and  the  analysis  of  the 
test  data,  the  following  conclusions  can  be  drawn. 

1.  The  total  sorption  of  carbon  dioxide  was  found  to 
be  greater  than  that  of  methane  and  the  sorption  of  methane 
in  turn  was  greater  than  that  of  nitrogen  for  all  samples. 

2.  The  ratio  of  the  total  amount  of  desorbed  carbon 
dioxide  to  methane  was  found  to  average  1.7:1,  but  this 
ratio  of  carbon  dioxide  to  methane  was  found  to  average 
3.4:1  at  the  end  of  the  first  27  seconds  of  the  desorption 
period. 

3.  No  regular  trend  in  total  sorption  capacity  was 
found  with  change  in  coal  rank. 

4.  Samples  from  Canmore  Mines  Limited  and  some  samples 
from  Crowsnest  Industries  Limited  showed  high  rates  of 
desorption  during  the  initial  period,  whereas  samples  from 
McIntyre  Mines  Limited  generally  showed  low  rates  of  desorp¬ 
tion  . 

5.  In  general,  the  quantity  of  gas  desorbed  Q  in  the 

first  minute  increased  linearly  with  /t  but  after  a  minute 

the  quantity  desorbed  can  be  represented  by  Q  =0^ ( l-Ae~X t) 

2 

where  A  =  tt  D .  This  expression  also  represents  the  increase 
a^ 

in  initial  rate  of  desorption  with  reduction  in  particle 
size.  Exception  to  this  was  noted  for  samples  from  a  highly 
fractured  and  sheared  seam,  in  which  case  there  was  little 
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difference  in  emission  rate  despite  differences  in  particle 
size . 

6.  High  values  of  the  index  were  observed, 

(indicating  seams  liable  to  outbursts).  In  general,  the  rela¬ 
tionship  between  the  AP  index  and  the  Qq_60  (methane)  was 
found  to  be  approximately  linear. 

7.  The  values  of  APq_^q  index  differed  from  layer  to 
layer  of  the  same  seam  (indicating  that  some  layers  are  more 
prone  to  outbursts  than  others) . 

8.  No  obvious  relationship  was  found  between  the  APrt 

0-6  0 

index  and  the  coal  rank  or  the  total  sorption  capacity. 


VOLUME  DESORBED  (CO2)  -  CU.  FT/TON  (50  pslg.  to  atmosphere) 
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FIGURE  5  VARITION  OF  TOTAL  SORPTION  CAPACITY  WITH  COAL  T>ANK 
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FIGURE  8  DESORPTION  OF  C02,  CH4,  N2  FROM  DIFFERENT  SIZES  OF  COAL 
Semple  2,  V/ILSON  SEAM,  Canmore 
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APPENDIX  I 


Calculation  for  desorbed  gas: 


V  =  volume  of  the  desorbed  gas  in  cc. 

A 

V2  =  total  constant  volume  of  the  apparatus  (cc.) 

=  1,780  cc . 

P1  =  atmospheric  pressure 

P  +  Ap  =  p2  =  atmospheric  pressure  +  pressure  due  to 

desorbed  gas 

AP  =  10  mv.  =  1.24  cm.  Hg  (by  calibration) 

(P1  +  AP)v2=  pi(v2  +  vx) 


Example  -  for  10  mv.  reading,  P1  =  70.0  cm 
71.24  x  1780  =  70.0(1780  +  V2 ) 

1780  x  1.24 


or 


V  = 
x 


70 


=  31.5  cc. 


Therefore  V  at  standard  atmospheric  pressure  =  31.5  x  -=-?■  = 
x  7  6 

29.0  cc.  The  conversion  factor  from  cc./gm.  to  cu.  ft. /ton  is 


106  1 
28.31  x  1.102 


32.1 


Therefore  total  volume  desorbed  by  a  6.5  gm.  sample  for  a 
10  mv.  signal  is 


x  32.1  =  144  cu.  ft. /ton 

Since  the  ordinates  of  Figures  6  to  16  are  to  a  scale  of  5  mv./ 
inch,  then  one  inch  on  the  ordinate  represents  a  yield  of  72 
cu.  ft. /ton,  measured  at  standard  atmospheric  pressure. 
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APPENDIX  II 

TRANSLATION  OF  DESCRIPTION  OF  AP  APPARATUS 
FROM  THE  LABORATORY  OF  CERCHAR 

Description  (see  sketch  no.  1  and  2) 

The  apparatus  is  composed  of: 

(1)  8  metal  cups  —  we  have  chosen  duraluminum  because 
of  its  low  density  —  intended  to  hold  8  samples  of  3  gm.  of 
coal.  Each  cup  is  provided  with  a  threaded  joint  which  assures 
the  seal  of  the  joint  with  the  glass  valve. 

(2)  8  glass  valves  with  a  flat  seat  at  the  bottom 
against  which  the  threaded  joint  of  the  corresponding  cup  is 
placed.  A  spherical  seat  grip  28/15  is  used  to  keep  each  cup 
and  valve  together.  By  turning  the  plug  of  the  valve,  the 
interior  volume  of  a  cup  can  be  connected  with  one  or  other 

of  the  two  manifolds  which  form  the  framework  of  the  apparatus. 

(3)  A  saturation  manifold  (front  manifold)  by  which 
the  samples  receive  either  helium  or  methane  through  the  purge 
valves  shown  on  the  right  side  of  the  sketch. 

(4)  An  exhaust  manifold  (rear  manifold)  joined  by  a 

line  to  a  two-stage  vaned  vacuum  pump  (the  exhaust  of  a  pump 

-2 

in  normal  service  is  limited  to  10  mm.  Hg)  and  also  con¬ 
nected  to  a  U-shaped  mercury  manometer  set  in  front  of  a 
millimetric  graduation  of  25  cm.  length. 

The  line  to  the  exhaust  pump  has  a  valve  for  the  entry 
of  air,  intended  to  facilitate  the  removal  of  the  cups  at  the 
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end  of  the  test. 

Each  gas  used  passes  through  a  dessicating  column  of 
colored  actigel  and  an  oil  splash  trap  which  allows  an  esti¬ 
mate  of  the  emission.  The  vaned  pump  is  protected  from  prod¬ 
ucts  which  might  contaminate  it  by  a  trap  of  liquid  nitrogen 
or  of  frozen  carbonic-acetone  according  to  their  availability. 
To  simplify  sketch  #1  we  have  not  shown  the  accessories  that 
have  just  been  mentioned,  which  are  shown  schematically  in 
sketch  #2. 

The  total  volume  offered  to  the  release  of  gas  after 

saturation  of  the  samples,  taken  to  the  mercury  in  the  zero 

position  of  the  manometer  (same  level  in  both  branches)  and 

deducting  the  actual  volume  of  the  coal,  is  equal  to  26  cm., 

3  3 

of  which  7  cm.  is  for  each  cup  and  19  cm.  for  the  exhaust 
manifold . 

Operating  Procedure 

3  grams  of  coal  to  be  tested  are  placed  in  each  cup. 

The  grain  sizes  adopted,  according  to  the  Russian  originators 

of  the  method,  are  as  follows: 

2-3  mm.  (-  -6+9  mesh  Tyler)  for  anthracites 
and  anthracitic  coals 

0.25  -  0.50  mm  (-32+60  mesh)  for  other  coals 
All  the  taps  are  turned  to  connect  the  cups  with  the 
exhaust  manifold,  and  the  pump  is  put  in  operation  for  1^ 
hours.  This  time  is  considered  sufficient  to  attain  an  equi¬ 
librium  of  degassif ication  under  the  chosen  operating  condi¬ 
tions.  The  taps  are  turned  to  an  intermediate  position,  the 
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front  manifold  is  exhausted  and  then  filled  with  helium  at 
atmospheric  pressure;  then  all  the  cups  are  connected  with 
this  manifold  and  the  operation  is  left  to  proceed  for  about 
10  minutes;  the  adsorption  of  this  gas  being  practically 
negligible  at  the  ambient  temperature,  a  longer  time  is  un¬ 
necessary. 

The  cups  are  then  isolated,  and  a  check  is  made  that 
the  manometer  is  in  proper  equilibrium  or  if  necessary  the 
corresponding  vacuum  is  re-established.  Each  cup  is  then 
connected  in  turn  with  the  rear  manifold  in  which  the  vacuum 
is  re-established  between  each  two  consecutive  operations. 

To  each  sample  corresponds  a  value  of  the  difference  in 
level  read  on  the  manometer,  expressed  in  mm.  Hg,  and  estab¬ 
lished  in  a  few  seconds. 

The  helium  is  removed  by  pumping  from  all  the  cups  at 
the  same  time  —  several  minutes  suffices  for  that  —  then 
methane  is  introduced  by  the  saturation  manifold,  under  at¬ 
mospheric  pressure  for  lk  hours,  at  the  end  of  which  time  it 
is  considered  that  equilibrium  of  adsorption  is  attained  for 
all  the  samples.  For  helium  as  for  methane  maintenance  of 
pressure  at  the  fixed  value  is  assured  by  a  lateral  connec¬ 
tion,  branched  ahead  of  the  purge  valve  and  passing  through 
a  splash  trap  common  to  the  two  gases. 

After  this  saturation  all  the  taps  are  closed,  the 
vacuum  in  the  exhaust  manifold  is  checked,  and  then  for  each 
sample  the  indices  P2  and  P^  are  measured  in  the  following 
way:  a  cup  to  be  tested  is  connected  for  10  seconds  with 
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the  exhaust  manifold,  the  tap  is  closed  and  the  difference 
in  level  P2  is  read  in  mm.  of  Hg;  at  the  end  of  35  seconds, 
i.e.,  45  seconds  after  the  beginning  of  the  test,  the  tap  is 
reopened  towards  the  exhaust  ramp  and  this  position  is  main¬ 
tained  15  seconds;  the  sample  studied  is  again  isolated  and 
the  total  difference  in  level  AP2  (P3?)  read  in  mm.  Hg,  cor¬ 
responding  to  a  desorption  period  of  60  seconds. 

Each  sample  is  characterized  by  the  speed  of  desorption 

indices 


AP0-10  P2  "  P1 


AP  =  P  -  P 

A *10-60  *3  *2 


AP0-60  AP0-10  +  AP10-60  P3  “  P1 


The  determination  of  P^  enables  knowledge  of  the  con¬ 
tribution  of  free  gas  which  is  non-adsorbed  on  the  sample  but 
simply  held  under  atmospheric  pressure  in  the  dead  volume 
formed  by  the  volume  of  the  coal  pores  and  the  intergranular 
spaces.  The  evolution  of  sorption  from  the  initial  moment  can 
thus  be  followed. 


Remark 

In  particular  for  samples  of  the  weakest  sizes,  the 
emission  of  coal  in  the  rear  manifold  is  avoided  at  the  begin¬ 
ning  of  exhausting  and  when  releasing  the  pressure  by  intro¬ 
ducing  a  small  pad  of  hydrophile  cotton  in  the  passageway  of 
each  of  the  taps. 
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Note 

A  letter  from  M.  de  Vergeron  of  September  16,  1965 

advises  that  "in  our  studies  concerning  the  characterization 

of  deposits  susceptible  to  sudden  outbursts  of  methane  we  take 

into  consideration  the  total  index  AP  =  AP.  , 

U  — 60 

Translated  T.  Patching,  Oct.  1965 
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APPENDIX  III  (a) 


Method  for  calculating  A  on  the  basis  of  Least  Squares  from 

i 

equation 

—  \  +* 

(5) 


Qi  =  Q.U  -  e'Ati) 


Supposing  =  f(A,t)  and  let  experimental  values  of  Q  be 


shown  as  Q . . 

1 


oo 

£ 


On  the  basis  of  Least  Squares  j_=1(Qt  -  Q^)  should  be  minimum 

t  [Q,  -  f(A,t.)]2 
i=l  1  1 

On  differentiating  with  respect  to  A  and  setting  this  equal  to 
Z0rO  [Qi-f(A,ti)]  f(A/ti)  =  0 


Qi  ax  f(x'ti)  _  i®1  fU(ti)  It  f(x-ti: 


00  00 

or  Q.  x  Q/„  t.e'Ati  =  Q/„  (i-e_Ati)  Q„  -t.  e_Ati 


00 

£ 


f ' (A)  =  - 


Qi 

t . 
i 

-  A  t . 
e  l 

-  Q 

00 

£ 

i=l 

Qi 

t.2  e- 
1 

At . 
l 

00 

£ 


00 

£ 


-At . v  -At 


-At.  .  -At 


-  t..(l  -  e  Xti)e~Ati] 


or  f '  ( A)  =  (-.f,  Q.  t.2  e  Xti)+Q  .f,  t . 2  e_Xti(l-2e  Xti) 

1  —  -L  11  oo  1 — 1  1 

Using  the  Newton  Rapson  formula  the  value  of  A  can  be  calculated 


>  -  *  f (A) 

l+l  Ai  FTT) 


as  follows: 
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An  A.P.L.  computer  program  was  written  to  solve  the  above 
equation  for  X,  and  was  given  the  name  "V  LAS QL AMD A  T" .  For 
this  program  the  following  symbol  changes  were  made: 

Q  +  V 

Q  +  VM 

00 

X  +  LAMDA 
t  -v  TIME 


Listing  of  function  V  LAS Q LAMDA  T 


VLASQLAMDAIU1V 

V  V  LASQLAMDA  T 

[  1  ]  U+FLD+i  +  /V*T*E)  -  (  VM*  +  /T*Ex  (  1  -  E<-* -LAMDA  *  T )  ) 

[  2  ]  n+FPLD-*-  (Q-<--  +  /7x  (T*  2  )  xE  )  +  □«-  VM*  +  /(T*2)xEx(l-2'*E) 

[3]  LAMDAU-LAMDA-FLDiFPLD 

[4]  ('LAMDA  ' \ LAMDA 1 ) 

L5]  +(1E~8>\ LAMDA1-LAMDA) /7 

[6]  +1,LAMDA+LAMDA1 

[7]  (» CORRECT  VALUE  OF  LAMDA  =  ' \LAMDA1) 

[8]  *  ’ 

[9]  ( ' TV E  VALUES  OF  V  USING  LEAST  SQUARE  APPROXIMATION 

VM*  (l-E<-*-LAMDAlxT)  ) 

V 
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APPENDIX  III  (b) 


To  calculate  values  of  A  and  A  from  equation 

Qt  =  Q„  (1  -  Ae~Xti) 


(6) 


Let  Q' .  =  f  ( A  ,  t . )  and  experimental  values  of  Q.  be  shown  as  Q. 

±  ±  t  1 

oo  o 

Then  [Qi  -  f(A,t^)]  should  be  minimum. 

By  differentiating  with  respect  to  A  and  A,  and  equating  to  zero 
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Q  •  —  = 
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Q.  M  = 
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<i>  i=i  QiQ^At.e'^i  =  Qjr-ST^ilQ^At. 
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-At. 
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(iii) 


oo 
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oo 

from  (ii)  il1  Qitie_xti  =  if1  t.^-*  (1-Ae'xti) 


(iv) 
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oo 

£ 


.  Q.e  Ati  =  .  e,  ( 1-Ae  Ati)  Q  e”Ati 

1=1  1  1=1  oo 


or 


or 


i=l  Qie'Xti  =  Q 


e  ^-At.  _  f 


A  = 


00 

e  At 
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1=1  1 


1  -  Q  A  .%  e‘2Xti 

oo  1=1 
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Now  substituting  value  of  A  in  (iii) 

...  -At.: 
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oo 
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i=l 
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-2At. 


i)  <i=i  Q«tie-Xti)  =  Q.  tie-Atii^1 
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e  -At 
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f(X)  =  (^e  2Xti)  Ati)-<2»iiitie  Ati)  <ilxe  2Xti> 


00  00  00  00 

+(iiltie  i) (Q«i=le_Xti) - (i=ltie_2Xti) (i=lQie_Xti) 


(i-itie"2Xti)I(i-iQiV"xti>-<Q-i-iV"Xti)] 


+  (i=itie  2Ati>  [(Q»i=ie"Xti)-(i=iQie  xti>l 


Differentiating  again  with  respect  to  A , 
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Using  the  Newton  Rapson  formula  the  value  of  A  can  be  calcu¬ 


lated  as  follows: 
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X 


i+1 


X  . 
1 


f  (X) 

FTT) 


A  computer  program  was  written  to  solve  the  above  equations 
for  X  and  A,  and  was  given  the  name  "V  LSSQ  T" . 


Listing  of  function  V  LSSQ  T 


NLSSQ [Q]V 
VV  LSSQ  T 

1 1  ]  S l«-(  (  +/TxE*2  )x  VM*+/E)  -  (+/27x#*2)x  +  /  Vx  E-*-*- LAMDA*  T 

[2]  FLD+( ( (+/E*2)x(+/VxT*E) )-( VM* + / Tx E ) x + /E* 2 ) +51 

[3]  S2«-(  2x  +  /(T*2)x£’*2)x(  (  Wfx  +/£)-(+/  yx£)  ) 

[4]  S3«-(+/2’x£*2)x(  (  7Mx  +  /2’xj?)-(+/7x2’x£’)  ) 

[  5]  S4+(+/£*2)x(  (  VM*+/ (T*2)xE)-(+/Vx(T*2)xE)  ) 

[6]  FPLD+-  ( 54+S3 ) -S2 

[7]  ('LAMDA  =  • iLAMDAl+LAMDA-FLD*FPLD) 

[8]  ->(1 E  8 >  |  LAMDA  1  -LAMDA  )  / 1  0 

[9]  ->1  ,LAMDA+LAMDA1 

[10]  (’raff  CORRECT  VALUE  OF  LAMDA  =  1 \LAMDA1) 
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